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Abstract 

We consider a system coupling the parabolic-parabolic Keller-Segel equations to the 
incompressible Navier- Stokes equations in spatial dimensions two and three. We estab- 
lish the local existence of regular solutions and present some blow-up criterions. For two 
dimensional Navier-Stokes-Keller-Segel equations, regular solutions constructed locally in 
time are, in reality, extended globally under some assumptions pertinent to experimental 
observation in [20] on the consumption rate and chemotactic sensitivity. We also show the 
existence of global weak solutions in spatially three dimensions with stronger restriction on 
the consumption rate and chemotactic sensitivity. 
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1 Introduction 

In this paper, we consider mathematical models describing the dynamics of oxygen diffusion 
and consumption, chemotaxis, and viscous incompressible fluids in M. d , with d = 2, 3. Bacteria 
or microorganisms often live in fluid, in which the biology of chemotaxis is intimately related 
to the surrounding physics. Such a model was proposed by Tuval et al. [20] to describe the 
dynamics of swimming bacteria, Bacillus subtilis. We consider the following equations in [20] 
and set Q T = (0, T] x R d with d = 2, 3: 

dtn + u • Vn — An = —V • (x(c)nVc), 

dtc + u-Vc- Ac = -k(c)n, in (x, t) G R d x (0, T], (1.1) 

k dtu + u ■ Vn — Am + Vp = — nV</>, V • u = 

where c(t, x) : Q T -)• R+, n(t, x) : Q T ->• M + , u(t, x) : Q T -)■ R d and p(t,x) : Q T ->• M denote 
the oxygen concentration, cell concentration, fluid velocity, and scalar pressure, respectively. 
The nonnegative function k(c) denotes the oxygen consumption rate, and the nonnegative 
function x{°) denotes chemotactic sensitivity. Initial data are given by (jiq(x), co(x), uq(x)). 
To describe the fluid motions, we use Boussinesq approximation to denote the effect due 
to heavy bacteria. The time-independent function <j> = cj>(x) denotes the potential function 
produced by different physical mechanisms, e.g., the gravitational force or centrifugal force. 
Thus, <j)(x) = axd is one example of gravity force, and <p(x) = (j>(\x\) is an example of centrifugal 
force. Experiments in [20] suggest that the functions fe(c) and x( c ) are constants at large c and 
rapidly approach zero below some critical c*. Hence, in [20], these functions are approximated 
by step functions, e.g., k(c) = k\9(c — c*) and x(c) = K20(c — c*) for some positive constants 
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K\ and K2- Also in [2], numerical simulation of plumes was obtained for the same species of 
bacteria in [20] in two dimensions. Furthermore, they assumed that the functions x( c ) an d 
k(c) are constant multiples of each other, i.e., x(c) = fj,k(c). 

The main goals of this paper are to show the local existence of smooth solutions in two and 
three dimensions with the general condition on the oxygen consumption rate and chemotactic 
sensitivity, and to demonstrate global existence of smooth solutions in two dimensions and 
weak solutions in three dimensions with appropriate assumptions of x{ c )i k(c), (f> and initial 
data. Here we mention the related works for the result in this paper. If we ignore the coupling 
of the fluids, we obtain the angiogenesis type system. The classical model to describe the 
motion of cells was suggested by Patlak[T7] and Keller-SegelpTJ [12]. It consists of a system of 
the dynamics of cell density n = n(t, x) and the concentration of chemical attractant substance 
c = c(t, x) and is given as 

nt = An — V • (n^Vc), 
act = Ac — Tc + n, (1.2) 
n(x, 0) = n (a:), c(x, 0) = c (x), 
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where x is the sensitivity and r" represents the activation length. The system in (|1.2p has 
been extensively studied by many authors (see [91 \TU[ [TBI [To] [21] and references therein). For 
the chemical consumption model by the cell or bacteria, we refer to the following chemotaxis 
model motivated by angiogenesis. 

n t = An — V • (nx(c)Vc), 

ct = -c m n, (1.3) 

n(x, 0) = no(x), c(x, 0) = cq(x). 

The global existence of the weak solution to the system in (|1.3p was obtained by Corrias, 
Perthame and Zaag[3j H] with a small data assumption of ||no|| d. The bacterial movement 
toward the concentration gradient model in the absence of the fluid, i.e., u = 0, was recently 
studied. When u = 0, x( c ) = X an d ^( c ) = c in ([l.ip . it was shown in [19] that there exists a 

uniquely global bounded solution if < x < — 7~, rn — n • 

6(d+ 1)||co||l«> 

If the flow of the fluid is slow, then the Navier-Stokes equations can be simplified to the Stokes 
equations. For the case x( c ) — Xi Lorz|14] showed the local existence of the solution for the 
Keller-Segel-Stokes system. In two dimensions, Duan, Lorz, and Markowich[6] showed the 
global existence of the weak solution to the Keller-Segel-Stokes equations with the small data 
assumptions on Co, <j) and the assumptions on the functions such that 

X(c) > 0, x'(c) > 0, fc'(c) > 0, ^ < 0. (1.4) 

In [13], Liu and Lorz showed the global existence of a weak solution to the two-dimensional 
Keller-Segel-Navier-Stokes equations with similar assumptions on k and x t° those in (|1.4p . 
The equation of n in (jl.ip could have been replaced by a porous medium equation, i.e., An is 
replaced by An m and the following Keller-Segel-Stokes system has been considered in [7]. 

dtn + u ■ Vn — An™ = —V • (x(c)nVc), 

d t c + u- Vc- Ac = -k(c)n, in (x, t) £ R d x (0, T}. (1.5) 

d t u - An + Vp = -nVcj), V ■ u = 0, 
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In [7], Francesco, Lorz and Markowich showed the global existence of the bounded solution to 
(jl.5p when m £ (|, 2]. In [13] . Liu and Lorz proved the global existence of the weak solution 
when m > | in three dimensions. For the Keller-Segel-Navier-Stokes system (II. ip . Duan, Lorz 
and Markowich[6] showed the global-in-time existence of the iif 3 (lR a! )-solution, near constant 
states, to (jl.ip . i.e., if the initial data ||(no — n^, cq, uq)\\ H 3 is sufficiently small, then there 
exists a unique global solution. 

As mentioned earlier, the aim of this paper is to obtain the local-in-time existence of the 
smooth solution in two and three dimensions and the global-in-time existence of the classical 
solution to in two dimensions under the minimal assumptions on the consumption rate 

and chemotactic sensitivity. Now we are ready to state our main results. The first result in 
this article is the local existence in time of the smooth solutions to (jl.lj) . Comparing with the 
result in [6], we show the local-in-time existence without smallness of the initial data. 

Theorem 1 (Local existence) Let m > 3 and d = 2,3. Assume that x> k £ C m (R + ) and 
k(0) = 0, HVVlli 00 < 00 f or 1 < Kl < m - There exists T > 0, the maximal time of existence, 
such that, if the initial data (n , c , it ) G H m ~ l (R d ) x H m (R d ) x H m (R d ), then there exists 
a unique classical solution (n, c, it) of (II. ip satisfying for any t < T 

(n,c,u) G L°°(0, t;H m ^(R d ) x FT n (IR d ) x H m {R d )), 

(Vn, Vc, Vu) G L 2 (0, t; i/" 1 ^ 1 ^) x i/ m (R d ) x H m (R d )). 

Remark 1 For simplicity, we denote 

\\(n(t),c(t),u(t))\\ Xm := \\n(t)\\ Hm -i {Rd) + ||c(t)|| H m(Hd) + \\u(t)\\ Hm{Rd) . 

We remark that if T is the maximal time of existence with T < oo in Theorem {J\ then 

]imsup\\(n(t),c(t),u(t))\\x m + [ II ("(*). <t), u(t))f x = oo, m > 3. 

Secondly, we obtain two blow-up criteria for the system (jl.ip depending on dimensions. 

Theorem 2 Suppose that x> k, <f> and the initial data (no,co,uo) satisfy all the assumptions 
presented in Theorem^ If T* , the maximal time existence in Theorem^ is finite, then one 
of the following is true in each case of two or three dimensions, respectively: 




Remark 2 Theorem^ can be interpreted as follows: If ||Vc||^oo < oo in two dimensions or 

io (I Ml Lp + ll^lli 00 )^ < oo in three dimensions, then the local solution can persist beyond 
time T, i.e., (n, c, it) G L°°(0, T + 5; X m ) n L 2 (0, T + 5;X m+1 ) for some 5 > 0. 
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The third main result is the global existence of the smooth solutions in the two-dimensional 
spatial domain M. 2 . Motivated by experiments in [20] and [2], we assume that the oxygen con- 
sumption rate k(c) and chemotactic sensitivity x( c ) satisfy the following conditions: 

(A) There exists a constant \x such that sup |x( c ) — ^( c )l < 6 f° r a sufficiently small e > 0. 

(B) x(c), k(c),x'(c),k'(c) are all non-negative, i.e., x(c), fc(c), x'(c), fc'(c) > 0. 

We remark that assumption (A) plays a crucial role in obtaining LlogL x ff 1 x L 2 type 
estimates. 

Theorem 3 (Global existence in two dimensions) Let d = 2. Suppose that Xi k, 4> and the 
initial data (no, cq, uq) satisfy all the assumptions presented in Theorem^ Assume further that 
X and k satisfy the assumptions (A) and (B) and <f> > 0. Then the unique regular solution 
(n, c, u) exists globally in time and satisfies for any T < oo 



Remark 3 If we approximate Heaviside functions using the smooth functions, then the exper- 
iments in jl20f satisfy the assumptions (A) and (B). Furthermore, the assumptions on <j) are 
satisfied by gravitational and centrifugal forces. Also we note that 2D numerical studies were 
performed under the assumption x(c) = fik(c) in Jl|/. 

Our final main theorem is on the global-in-time existence of weak solutions in three dimen- 
sions. The notion of a weak solution of (jl.lj) is detailed in section 4 (see Definition [5]) . For 
existence of global weak solution, we need similar restrictions on k(c) and x{ c ) as m Theorem 
El More precisely, compared to (A), we impose a slightly stronger assumption, denoted by 
(AA), which is given as follows: 

(AA) There exists a constant [i such that x( c ) ~ ^k(c) = 0. 

We are ready to state our last main result. 
Theorem 4 Let d = 3 and (no,co,uo) satisfy 



Assume further that x an d k £ C 1 (]R + ) satisfy the assumption (AA), (B) and k(0) = and 4> 
satisfies </> > and ||V'0||l°° < oo for < \l\ < 2. Then a weak solution (n,c,u) exists globally 
in time. 

The rest of this paper is organized as follows. In Section 2, we prove local-in-time existence 
of the smooth solution for the two and three dimensional chemotaxis system with incompress- 
ible Navier-Stokes equations and obtain some blow-up criteria for the solution. In Section 3, 
we show the global in time existence of the regular solution in two dimensions. In Section 4, 
we establish the existence of a weak solution in three dimensions. 



(n, c, u) G L°°(0, T; H m ~ 1 (R 2 ) x H m (R 2 ) x H m (R 2 )), 
(Vn, Vc, Vu) G L 2 (0, T; H m ~ 1 (R 2 ) x H m (R 2 ) x H m (R 2 )). 



n G L 1 (R 3 ), c G H 1 ^ 3 ) n L°°(IR 3 ), u G H(R 3 ) 





(1.8) 
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2 Local existence and blow-up criterion 



2.1 Local existence 

We first consider the chemotaxis system coupled with the Navier-Stokes equations in two and 
three dimensions. We show the local existence of solutions (n, c, u) in i? m_1 x H m x H m space 
with m > 3. 

Proof of Theorem 1. We construct the solution sequence (n J , c 7 , u J )j>o by iteratively 
solving the Cauchy problems on the following linear equations 

( d t ni +1 + v? • VnJ +1 = An 3+1 - V • (x(c j )n 3+1 Vc 3 ), 

d t ci +1 + v? ■ Vci +1 = AcJ +1 - k(d>)n^, (2.1) 
k d t u j+1 + u j ■ Vm j+1 + Vp> +1 = Au j+1 - n j V(j), V • u j+1 = 0. 

We first set (n°(x,t), c°(x,t), u°(x,t)) = (no(x), co(x), uo(x)). Then, using the same ini- 
tial data to solve the linear Stokes type equations and the linear parabolic equations, we obtain 
u l (x,t), c (x, t) and n 1 (x,t), respectively. Similarly, we define (n J (x, t), c 3 (x,t), u^(x,t)) iter- 
atively. For this, we presume that c? and n J are nonnegative and show the existence and the 
convergence of solutions in the adequate function spaces. We show the nonnegativity of c J and 
■n? at the end of the proof. 

To prove the conclusion, i.e., to obtain contraction in adequate function spaces, we show the 
uniform boundness of the sequence of functions under our construction via energy estimates. 
• (Uniform boundedness) We here show that the iterative sequences (n J , c 3 ,u J ) are in X m := 
H m ~ l x H m x H m space for all j > 0. Observing that 

\\d a (u j n j+1 )\\ L 2 < C{\\u j \\ La o\\n j+1 \\ Hm -i + ||n J ' +1 || Lco \\u j \\ H m-i), 

\a\<m—l 

11x^)^+^11^^ <C\\nJ +1 \\ Hm _ L (l + H^ll^) , 
||fc(o ? ')?i J |L m _ 1 < ClMI x (l + llc^ll™ M , 



we have the following energy estimates: 
(i) The estimate of n J+1 



Ii*. || n i+i|| 2 + ||v^' +1 || 2 <C\\u j \\ \\n j+1 \\ IIVn j+1 ll 

2 dt uH" 1 ^ 1 II Wh™ 1 ^ 1 — II \\l°° II II H" 1 ^ 1 II \\H m ~ 



+C '\\jjm-i || n "' +1 ||jjm-i ||Vn- ?+1 ||jj m _i + C ( 1 + || cJ ||jjm) || rj, ' ? ' +1 llffm-i ||Vn : ' 



7-4-1 1 1 

|ffm-l 



< C (l + \\v?f Hm + \\^ m ) \\n^% m ^ + \ HVn^H^ . (2.2) 



(ii) The estimate of c J+1 
1 d 



^'-i- 1 !! 4- ll\7^' +1 ll <n\\iii\\ \\ri +1 \\ ll\7^' +1 ll 

2 dt -H" m II II H m — 1 1 1 1 L 00 1 1 II H m II II H 71 



A-C WijW ll\7/^ +1 ll A- C ( 1 -I- 1 1 r^' 1 1 m I IItj^'II IIV^ + 1 

~ ii \\H m II ii/f m II \\H m v II \\H m i II iih" i_i ii uj?" 
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< C (l + ll^'H^ 1} ) ll^'H^ + C \\v?f Hm \\c J+1 \\ 2 Hm + \ ||Vc? +1 ||^ m . (2.3) 
(iii) The estimate of u^ +1 

\jt H nJ+1 fe m + ^ c ll^'ILoc lk' +1 ll^ 

+C7||u J '|| flW ||V^ +1 || LOO ll^l^+Cll^H^ ||V^ +1 ||^ m 

< C|K'fc m ||^' +1 t m + C|H|^ m _ 1 + i ||V^ +1 ||^ m , (2.4) 

where standard commutator estimates are used. We show that there exists a constant M > 
such that, for any j, the following inequality holds for a small time interval [0, T] (T will be 
specified later): 

sup ( ! + lie 7 ')! „ m + 11^ II 

0<t<T V 

\jjm— 1 "T" || v ^ ||//m ~r || v ix ||j^ 

Here M is a number with M > 4(||ra ||#m-i + ||co||# m + ||^o||#m). 

We prove (|2.5p via an inductive argument. Suppose (|2.5p hold for j < i. If we add (|2.2p . 
(|2.3p . and (|2.4p and use Young's inequality, then we have 



+ / IIV^'ll* + ||V^||' + IIV^'IlL, dt < M. (2.5) 



o 



d ,\\ ,'4-1 II 2 . || ,'4-1 l|2 . || ,--4-1 1 1 2 N . ||„ ,'4_1 ||2 ||„ ,-4_l ||2 . ||„ ,'4-1 i|2 



(\\n i+1 \\ + llc m H +llti i+1 H ) + \\X7n i+1 \\ + IIVc i+1 H + IIVu i+ H 
Ml II H 777 ^ 1 II I Iff™ II \\H m ' II UH 777 ^ 1 II \\H 777 II uH" 



^ C ^1 + II^H/fm + || c *||/fm) 1 1 71-4+1 1 1 H m ~ 1 + C (^1 + ll C *ll/i" m ) ll^lliY™- 1 
+C ||u ||^- m ||c || jj m + C \\ u \\ flm 11"^' \\ H m ^ \\ n \\ H 777 - 1 

<C(1 + M + M m ) lln^ 1 !!^.! + CM \\c i+1 \\ 2 Hm 

+CM ||ti i+1 ||^ m + C(l + M m - l )M + CM. 
In the last inequality, we use the induction hypothesis. Hence, we get 

j t i\\n l+1 \\ 2 Hm -, + V + % m + h l+l f Hm ) + ||Vn ?+1 ||^ m _ 1 + ||Vc l+1 ||^ m + ||Vn i+1 ||^ 



< C(l + M + M m )(||n i+1 ||^ m _ 1 + ||c i+1 ||^ m + ||n i+1 ||^ m ) + C(l + M m ^)M. (2.6) 

We choose time T such that max{C(l + M + M m )T, C(l + M m_1 )T} < \. Then from 
Gronwall's inequality, we have (|2.5p . In short, we have (p? + , d> + , ii J+1 ) G L°°(0, T;X m ) and 
(VnJ +1 , Vc-' +1 , Vii J+1 ) G L 2 (0, T;X m ) and the uniform bound ([22]) for small T > 0. 

Also if we multiply {c^ +l ) q ~ 1 on the both sides of the second equation of f)2. 1 1) and integrate 
over spatial varaibles, then we obtain 

^II^IIL + ^llv^fiil^o. 

Thus, the L°° norm of cP +1 is uniformly bounded, which implies that x( c ' 7 ) an d k(c*) are uni- 
formly bounded for all j > 0. 
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• (Contraction) The estimate of this part is similar to that of the previous one. For con- 
venience, we denote 8fi +l := — P ■ Subtracting the j-th equations from the (j + l)-th 
equations, we have the following equations for <5n J , <5c ?+1 and 5u^ +1 : 

( d t 5ni +1 + v? ■ V5ni +1 - A5n^ +1 = -Su j ■ VnP - V • (x(c>)5ni +1 Vc>) 

-V • (x(c>)n : >Vc i ) + V • (x(c 7 '" 1 )n , 'Vc'- 1 ), 

d t 5ci +1 + u> ■ V5ci +1 - Add> +1 = -5u^ ■ Vc^' - k(c*)6n> + (fc(c*) - k(c^ l ))n^\ (2.7) 

dt8u> +1 + v? ■ V5ui +1 + V8p> +1 - A<)V +1 = -Sv? ■ Vui - 5n^Vcp, 
k V • 5u j+1 = 0. 

(i) The estimate of <W +1 . 

Using the following standard commutator estimates 

Y J [d a {v? ■ V6n j+1 )d a 5ni +1 - (u j ■ Vd a 5n j+1 )d a Sn j+1 ] 

\a\<m—l 

< C(||V^|| L oo|| ( 5nJ +1 ||^ m _ 1 + ||u i || i y m -i||V^' +1 || L o=||^ +1 || jFfro - 1 ), 
we have the following estimate: 

\j t W 5n3+1 \\H^ + l|V<Sn i+1 |&m-i < C||V^|| LO o|| ( 5^ +1 ||^ m _ 1 

+C\\u j \\ Hm -i\\6n j+1 \\ Hm \\5n j+1 \\ Hm -i + C\\5u j \\ L o \\n j \\ Hm -i\\V6n j+1 \\ Hm -i 
+C\\8u j \\ Hm -i\\n j \\ L ~\\V8n j+1 \\ Hm -i + 0(11^11^ + ||c J '||^ 1 )||<Jn J ' +1 || Hm -i||V<yn J ' +1 || Hm -i 
+C\\n j \\ Hm -i(l + ll^'H™- 1 + || C ? - 1 ||^- 1 )|| ( 5 C J '||^||V^ +1 ||^-i. 

< C( M'ZxJ'" H^^. ^Il-zx-^ M^r^ ) ll^^"- 1 - 1 ^ H-C||r5^" 1 1 1 1 1 1 x -K C ( 1 1 c^" 1 1 ^ 1 1 1 1 1 > ) 1 1 ^ 1 1 1 x 

+£711^111^(1 + + p- 1 !!^)!!^!!!™ + i||v^+ 1 ||| w _ 1 . 

(ii) The estimate of 6c? + . 

l±\\6c? +x \\ 2 Hm + ||V^ +1 ||| ra < C||^|| i o ||^+ 1 || ifm ||V^'+ 1 ||^ 

+C\\u j \\ Hm -i\\6c> +1 \\ H m || V6d> +1 \\ H m + C\\5u j \\ L °o \\c> \\ H ™\\V 5c> +1 \\ H ™ 
+C7||c?'|| jffm ||^"||^- 1 ||Wc 7 ^ 1 || jF r- H- C(l + ||c^' l^" 1 ) ||^-! || Wc 7 ^ 1 1| ^ 
+C\\n^ 1 \\ Hm - 1 (l + H^ll^- 1 + ||^- 1 ||™- 1 )||5^||^||V5^ +1 ||^ 

< C||^||lr-ll^- , - :L |l^- CII^H^II^- 7 "!!^- C(l -r- ||^||^ _:L) )||^^||^ m _a 

+CK'- 1 11^(1 + \W\\%Z~ 1] + {tf-H^WJfir* + ^l|V^ +1 |||, m , 

where we used the Mean Value Theorem for the last term. 
(Hi) The estimate of 5u^ +1 . 

^¥u j+1 f Hm + ||V5u J+1 ||^ m < C||V^|| L -||(5n J ' +1 ||^ m 
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+ C\\ U j \\ H m || V^n^ 1 || £ oo ||^' +1 \\ H m + C\\5u j || £ oo || W ||tfm-l ||5u J ' + 1 ||#m+l 
+ C\\5u j \\ H m-l\\Vu j \\ L oc\\5u j + 1 \\ Hm + l + CllfWll^m-lHfo^ll^m+l 

< C||^V™ll^ j+1 Hlr™ + C\\ui\\ 2 Hm \\5vPf Hm +C\\Sn j \\ 2 Hm ^ + ^ ||<Su J ' +1 ||| m+1 , 

where similar standard commutator estimates are used as in the case of <W +1 . Using Young's 
inequality, we have 

^(ll^ +1 |||, m - 1 + ||<5c J ' +1 ||^ m + ||<5^ +1 ||^ m ) + ||V^ +1 ||^ m _ 1 + ||V^ +1 ||^ m + ||V^' +1 |||, m 

< C\\SnP +1 \\%^-^ +C\\5c> +1 \\ 2 Hm + C||<5^' +1 ||^ m +C\\5ni\\ 2 Hm - 1 +C\\5d\\ 2 Hm +C\\5u ] \\ 2 Hm , 

where C depend on the H m ~ l x H m x H m norm of (nP, c? , vP) and (nP~ , cP , and the 

maximum values of an d . Gronwall's inequality gives us 

sup (||W +1 |||, m _ 1 + ||^ +1 ||^ m + \\5ui +1 \\ 2 Hm ) 

0<t<T 

<CTexp(CT) sup \\ 2 Hm ^ + \\8c? \\ 2 Hm + \\5u j \\ 2 Hm ). 

0<t<T 

From the above inequality, we find that (n 3 , e 7 , u J ) is a Cauchy sequence in the Banach space 
C(0, T; X m ) for some small T > 0, and thus we have the limit in the same space. 

• (Uniqueness) To show the uniqueness of the above local-in-time solution, we assume that 
there exist two local-in-time solutions (ci(x,t), ni(x,t), U\(x,t)) and (c 2 (a;, t), ri2(x,t), u 2 {x,t)) 
of (|l.ip with the same initial data over the time interval [0, T], where T is any time before the 
maximal time of existence. Let c(x,t) := c\{x,t) — c 2 (x,t), n(x,t) := n\(x,t) — n 2 (x,t), and 
u(x, t) := u±(x, t) — U2(x,t). Then (c, n, u) solves 

' dth + u\ ■ Vn — An = —u ■ Vn 2 — V • ((%(ci) — x{ c 2))niVc{) 
-V • ( X (c 2 )nV Cl ) - V • (x(c 2 )n 2 Vc), 

d t c + ui • Vc - Ac = -u ■ Vc 2 - (fc(ci) - k(c 2 ))ni - k(c 2 )n, in (a;, t) £ R d x (0, T). 

dtu + u\ ■ Vtt — An + Vp = —u • Vu 2 — nVcfi, 
I V ■ u = 0, 

(2.8) 

Multiplying n to both sides of the first equation of (|2.8p and integrating over M. d , we have 

T^INl! 2 + II v ^IIl2 ^ ||"Sn2||L2||Vn|| L 2 + \\(x(ci) - x(c2))niVci||£2||Vn|| L 2 

+ ||x(c2)nVci|| L 2||Vn|| i2 + ||x(c2)n 2 Vc|| i 2 1| Vn|| L2 
< C||u||£ 2 Hr^llf/x) + C||c||^ 2 ||ni|||oo ||Vci H^oo + C||n||^ 2 ||Vci||^oo 
+C||Vc||^ 2 ||n 2 |||o +e||Vn||| 2 . 
Multiplying c and — Ac to both sides of the second equation of (j2.8|) . we have 

^^ll^lli 2 + H^lli 2 < ||«Q2||i,2||Vc||x,2 + C||c||| 2 ||ni|| L oo +C||n||x,2||c|| i 2 
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1111O1 o 1111O1 1111O 1111O 11O 

< C||m|| l2 ||c2|| i ,oo + C||c|| L2 ||ni||icx) + C||n|| L2 + C||c|| L2 + e||Vc|| L2 , 

and 

-^||Vc||| 2 + ||Ac||| a < \\ui ■ Vc|| i2 ||Ac|| L2 + \\u • Vc 2 || i 2||Ac|| L2 
+ ||A;(ci) - /c(c2)|| L 2||ni|| L oo||Ac|| L2 + C||n|| L2 1| Ac\\ L 2 
< C||V5||| 2 ||m|||oo +C||u|||a||Vp2||x ( oo +C||c||| 2 ||ni||| 00 +C||n|||a +e||Ac||| 2 . 
Multiplying u to both sides of the third equation of (j2.8|) and integrating over we have 

ld_ 

2dt' 

< C||{t||^ 2 ||« 2 |||oo +C\\n\\ 2 L2 +C\\uf L2 + e||Vu||| 2 . 
Summing the above estimates, we obtain 

1 d 

2df ' 



n~INli 2 + ll^^lli 2 — II^IIl 2 II^IU 00 1| Vn||^2 + C||n||£2 \\u\\ L 2 



;^(l|n||| 2 + ||c||| 2 + ||Vc||| 2 + ||«||| 2 ; 



< C(|| I Vci| + |Vc 2 | |||oo + ||ni + «2||l°° + ll n illi°° ll^cilli 00 + ll c 2|||°° + II l^i| + 1^1 |||°° + 1) 

x(ll^lli 2 + l|c|li 2 + ||Vc||2 2 + ||n||2 2 ). 

Since all L°° norms of (n,,Cj,Uj) are controlled by H m x H m x H m norm of (nj,Cj,iij) with 
m > 3, and the initial data of (c, n, u) are all zero, (c, n, u) are all zero for T > 0. That implies 
the uniqueness of the local classical solution. 

• (Nonnegativity) For completeness, we briefly show that n 3 and c 1 are nonnegative for all 
j. To use induction, we assume c 3 and v? are nonnegative. If we apply the maximum principle 
to the equation of c? +l in (12. ip . we find that c ?+1 is nonnegative (/c(c ? )n J is nonnegative). Let 



_ n i+i n j+i < o 
n j+1 > 0. 



= n 


3+1 
+ 


3+1 

— n_ . 


where 






hi n i4 


1 > 


H 


1° 




^0, 



j+l I ^ u j+1 

n, = < . , , n_ 



Recall that the weak derivative of n{_ +1 is — Vn J+1 if n{_ +1 < 0, otherwise zero. It holds that 



o 



d t ni + \ni +1 )„dxds = i(||(n'' +1 )_(0)IU 2 - ||(^ +1 )_(0)|| L2 ) 



\ jf ||(^ +1 )_||i 2 + ||V(^ +1 )_||i 2 d S < C jf* ||(^ +1 )_||i 2 ||Vc'||l 00 + l||V(^ +1 )_||f 2 d S . 



since rrL +1 , (9 t n J+1 G L 2 (0, T; L 2 (]R 2 )) (see e.g. p2]). Now multiplying the negative part (n J,+1 )_ 
on both sides of the first equation of (|2.1|) and integrating over [0,i] x M 2 , we have 

Using Gronwall's inequality, we have 

||(^' +1 )-(t)||| 2 < ll(n J+1 )-(0)|l! 2 exp (c£ IIV^'II 2 .^ . 

Since the initial data is nonnegative, we conclude that n- 7 " 1 " 1 is nonnegative. This completes 
the proof. rj 
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2.2 Blow-up criterion 

Next, we observe a blow-up criterion for the fluid chemotaxis equations. 

Proposition 1 (A Blow-up criterion) Suppose that \, k, <f> and the initial data (no,co,uo) 
satisfy all the assumptions presented in Theorem^ If T < oo is the maximal time of existence, 
then 

r T . . 

/ /.Ml II—, / ,\ ||9 \ .. 

oo . 



\\Vu(t)\\ Loa{wd) + \\Vc(t)\\ 2 Loom )dt 

Proof. At first, we consider the L 2 estimate of n. Multiplying n to both sides of the equation 
of n and integrating, we have 

^IMlia + ll Vn lli2 < C\\x(c)nVc\\ L 2\\Vn\\ L 2. 
Since x is continuous and c is uniformly bounded until the maximal time of existence, we have 

C||x(c)riVc|| L 2||Vn|| L 2 < |||Vn||| 2 + C||Vc|||oo||n|||a. 
For the estimates of c, we use the calculus inequality 

||V(u • Vc) - (it • V)Vc|| L 2 < C\\Vu\\ L o°\\Vc\\ L 2. 
Multiplying —Ac to both sides of the equation of c and integrating, we obtain 

— llVcUla + ||Ac||| 2 < C||Vn||^||Vc||| 2 + C\\(k(c)n)f L2 + \\\Acf L2 . 

For the equations of u, multiplying —Au to both sides of the equations and integrating by 
parts, we have 



1 d 
~2dt 



|Vn||| 2 + || Au||| 2 < C || Vw|| L oo || Vm|| L 2 + C ||n|| L2 || Aii|| L2 • 



Collecting all the estimates, we obtain 

j t (iHlls + ||Vc||2 2 + ||Vu||| 2 ) + {\\Vnf L2 + ||Ac||£ 2 + ||A«||| 2 ) 

< C (||Vc||ico + ||Vu||£oo) (||n||| 2 + ||Vc||| 2 + ||Vu||£a) . 
From Gronwall's inequality, we have 

sup (\\n\\ 2 L2 + ||Vc||2 2 + ||Vn||2 2 ) + f (||Vn||| 2 + ||Ac||^ 2 + || A«||| 2 ) dt 

Jo 

< C(\\n \\ 2 L2 + ||V Co ||| 2 + ||V«o||i a )exp (£ ||V«|U=c + \\Vc\\ 2 Loo dt\ . 

Note that |M|l°o(o,t ; l 2 ) an d \\^ n \\L 2 (o,T;L 2 ) are uniformly bounded if J Q T || Vit||x,°° + HVcH^oodt 
is bounded. Moreover, we see that n £ L%L^° and Vn 9//2 £ L 2 L 2 for all 2 < q < oo. Indeed, 



dt M ^ + 



_ 1 
Vn 2 


2 W 


InVcVn 9 " 1 ! 






2 



Vn2 



2 

L 2 
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From the above inequality, we have ||n(i)||xg < C, where C is independent of q. Letting 
q — > oo, we have n G L^Lf . 

Next, we consider the estimate in the space (n, c) G -ff 1 x if 2 . We have 

^ HVnllia + HAnHla < C \\Vu\\ LOO ||Vn|| 2 2 + C ||Vn|| L2 ||Vc|| LOO ||V 2 n|| i2 

+C ||n|| LOO || Ac|| L2 ||V 2 n|| L2 + C \\n\\ LOO ||Vc|| £00 ||Vc|| L2 ||V 2 n|| L2 . 
From Young's inequality and Gronwall's inequality, we have 

/•T 

sup||Vn|| 2 2 + / ||V 2 n|L 2 ^ 

JO 

< ^llVnoll^a +C||n|| ioo(0)T . z ,oo) ^ l|Ac|| 2 2 dt+ ||Vc|| ioo(0iT . L2) ^ HVcH^dt 

xexp^y ||Vu||x,oo + ||Vc||x,oodt 
Hence, n G H^.L^° n H^L 2 . For the -ff 2 estimate of c, we have 

~ ||Ac|| 2 2 + ||VAc|| 2 2 < C \\Vu\\ LOO ||Ac|| 2 2 + C ||Au|| L2 ||c|| LOO ||VAc|| i2 

+C ||Vc|| L 2 \\n\\ LOO ||VAc|| L 2 + C \\Vn\\ L 2 ||VAc|| L2 . 

By Gronwall's inequality, we have c G H%Lf n H 3 L \. Similarly, u G H%L n if 3 ^ 2 . Then, 
we consider the estimate in the space (n,c, u) G if 2 x if 3 x H 3 . Proceeding similarly to the 
above, we obtain 

^^ll n ll^ 2 + ll Vn ll// 2 ^ ^IMU°°ll ra llir 2 ll^ n ll.ff a + C||Vtt||£,oo||n|| H i||Vn|| H 2 

+ ^l|Vn||| 2 +C|| X (c)nVc||^. 
In the above, the last term can be controlled by 

\\ X (c)nVc\\ H 2<C\\n\\ H2 \\ X (c)Vc\\ H2 , 

and 

||V 2 (x(c)Vc)|| L2 < C ||V 3 c|| L2 + C ||V 2 c|| L2 \\Vc\\ Loc + C ||Vc||| 6 . 

We already obtained c G H^L^° n H 3 L 2 . Hence, if we use Young's inequality and Gronwall's 
inequality, we have 

sup 1^11^2 + J ||Vn||^2 dt < ||no||# 2 exp + J ||Vu||l°° + HVcH^cZ^ . 
Similarly, we estimate c as 

kills + ||Vc||^ 3 < a||V«|| L cx.|| c ||^3||Vc||^3 + C\\u\\ H 4Vc\\ L o* ||Vc|| h3 



-— llrl' 2 
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+C\\(k(c)n)f H2 + \\\Vcf H3 . 

We can control the term ||(&(c)ra)||^ 2 °y C(||c||^ 2 IMIif 2 + ll c ll#i II^ c IIl°° IMIh 2 )- For the 
estimate of u, we have 



T^IMI^a + ll Vn llH 3 - C|l Vu l|L°°IMIi/ 3 ll Vn ll#3 + ^ 



Using Gronwall's inequality, we have (c, u) £ (H%L%° n i^L?) x (H^L™ n H^L 2 ). Let us 
consider if m_1 x iJ m x estimates. The case m = 2, 3 and 4 are proved in the above, hence 
we consider the m > 5 case. Taking 9" (|a| < m — 1) and multiplying <9 a ra to both sides of 
the equation n and integrating and summing, we have 

_■ (ML 

+-||V?i||^ m -i +C\\x(c)nVc\\ 2 Hm - 1 . 

We already obtained the estimate for the case m = 4, thus HVcII^oo^TiL 00 ) is bounded. Hence, 
we have 

||Vx(c)||flm- a < C(l + ||Vc|| L =o)||Vx'(c)||flm-3. 

Using the classical product lemma on each step of iteration, we can control 

\\V X (c)\\ Hm -2 <C(l + \\Vc\\ Loa ) m - 1 . 

Then we have 

||x(c)nVc|| Hm _i < (7(1 + ||c||irm + ||Vc||£») \\n\\ Hm -i 
using the product lemma. For the H m estimate of c, we proceed similarly to have 

1 d 
~2dt 



Lao ||c||#m || Vc|| JJm + C||u||/fm || Vc||ioo ||Vc||j^r, 



+C||(fc(c)n)||l P n-i+i||Vc|||p n . 

As is shown for the term ||x(c)nVc||^ m -i , we control the term ||(A;(c)n)||^ m _i by (7(1 + 
||Vc||^T 2 ) ||n||^ m -i. For the estimate of u, we have 



^^IMlIf™ + ll^lllf" 1 - C|| VTi||ioo ||it||jym ||Vn|| /J - m + i 



Thus, by collecting all the above estimates and using Gronwall's inequality, we have (re, c, u) £ 
(H™~ l LfC\H™L 2 t ) x (H™L™nH™ +1 L 2 ) x [H^Lf n H™ +1 L 2 ). This completes the proof. n 

We are ready to present the proof of Theorem [2j 
Proof of Theorem 2. In the proof of Proposition 1, we notice that ||Vc||l°o is solely 
responsible for re £ L^Lf 3 and Vn £ L 2 L 2 . Indeed, 



d „ 
dt 11 1 



£3 + || Vn||Ja < C f |nVcVre| dx < C \\Vc\\ 2 Lao \\n\\ 2 L2 + ]- \\Vn\\ 2 L 2 . (2.9) 
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This implies u G L^Lf and Vu G L^Ljf by 



2r2 



^ II l|2 , ||V7 1 1 2 ^ r,\\ II || || 

-j t IMIl 2 + II v «IIl2 < C||n|| L 2||n|| i2 . 
Moreover, we have n G L%L^° and Vn'/ 2 G £ 2 L 2 for all 2 < g < oo; 



(2.10) 



NIL* + 



+ 






nVcVn^ 1 






l 2 h 





Vn' 



2 

L2 



Next, we see that Vc G L 2 x Lf and V 2 c G £ 2 £ 2 . Indeed, 

^ l|Vc||^a + ||V 2 c||^ 2 < C||Vc|| Loo 
We first consider the two-dimensional case. 



u \\l 2 IK C \\l 2 ^ IHI-L 2 IK C \\l 2 ' 



(2.11) 



• (2D case) For convenience, we denote vorticity as w := V x u; that is w = d\U2 — c^l 
in two dimensions. Next, we consider the vorticity equation 

ujt — Aw + iiVw = — V riVcf), 

where V^n = (— d^n, din). We note that w G L^Lf and Vw G L 2 L 2 , since 



— Hwll^a + ||Vo;||^2 < C||Vn|| L 2 |M| i2 



| L 2 T || VW|| L ! 

'2 roo „^ V72, 



(2.12) 



Furthermore, we observe that Vw G L^L^ and V^w G L 2 L 2 . Indeed, testing —Aw, we get 
d „_ ,,2 



,|Vw|| f 2 + ||V 2 w||f 2 < llitl 
dt 1 L 11 llL ~ 



L 4 ||Vw|| L 4 ||Aw|| L 2 + ||Vn|| i2 ||Aw|| i2 



<C\\u\\ 2 L2 \\Vu\\ 2 2 \\Vu\\ 2 2 \\V 2 uj\\ 2 r2 + || Vn 



L 2 



L 2 



\L 2 



\L 2 



I Awl 



Therefore, via embedding, we have 

r, , /" r „ , / ,T „ 

/ ||Vn|| ioo <ii < / ||Vu||^2 < C / ||w||^2 eft < oo. 
Jo io io 

This completes the proof of the 2D case. 

• (3D case) We will show this case by contradictory arguments. We suppose that the con- 
dition (|1.7p is not true. We first recall the vorticity equation 

wt — Aw + uVw = wVu — V^nV0. 

Under the condition (|1.7p we have w G L 2 Lf D and Vw G -£/ 2 £ 2 as follows. We denote Q* = 
M 3 x (T* - <$,*) for T* - 5 < t < T*. For any given p, q satisfying 3/p + 2/g = 1, 3 < p < oo, 
we choose l,m such that l/p+ l/l = 1/2 and l/q + 1/m = 1/2. We then remind that, due to 
the Gargliardo-Nirenberg's inequality, 



it ,i 



. in 

< C \\u\\ r 2 IIVull 1 ^ , 2 < / < 6, - + 



I m 
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where 6 = (6 - l)/2l and 1-6 = (31- 6)/2Z. Then we have 

^ IMl£a + HVwH^a < \\u\\ LP \\u\\ L i \\Vuj\\ L 2 + C \\n\\ L 2 \\Vu\\ L 2 
Next, integrating in time over (T* — S,t), 



\\u(t)\\h + / < M T * ~ + C IMIl™ IMI l ^ II Vw|| L 2, 2 + C ||n|| L2 , 2 || V W || £2 ,2 

< ||w(T* - 5)11^2 + IMIi£9 ||w|| L 2,oc 1 1 Vw || + C ||n|| L 2,2 HVwIl^a . 
Note that 9 > 0. By Young's inequality, we have 
IK*)llis + / l|Vw||^ < \\u(T* - S)\\ 2 L2 + C \\u\\%, q ||w||? a ,cc + C \\n\\ 2 T 2,2 + J ||Vw|| 2 2 , 2 . 

Jf* —8 x,t x,t x,t 2 x,t 

Since ||tt]|^9 can be sufficiently small in (T* — 5, T*) x R 3 by decreasing 5, we have 

IK*)lli2 + \f ||Vu;||£ 2 < ||u,(T* - 5)\\ 2 L 2 + C ||n|| 2 2 , 2 , (2.13) 

Z J T *_S *,t 

which is bounded by (|2.9p . Since f is arbitrary for all t < T*, this estimate is uniform. 
Next, we observe that V 2 c G L^Lf 3 and V 3 c G L 2 L 2 . Indeed, we estimate 

4 ||V 2 c||* 2 + ||V 3 c||^ 2 < C(||Vn|| L2 + ||Vc|| Loo ||Vn|| L2 ) ||V 3 c|| L2 + ||u|| L «, ||V 2 c|| L3 ||V 3 c|| L2 

< C(||Vn|| L2 + \\Vc\\ LOO \\Vu\\ L2 ) ||V 3 c|| L2 + ||w|| L2 ||V 2 c||| 2 ||V 3 c||J 2 , 
and use (|2.9|) . (|2.10|) . Similarly, we show that n G Lf>H]. n £ 2 # 2 by estimating 

^ ||Vn|| 2 2 + ||V 2 n||' 2 < ||u|| £6 ||Vn|| L3 ||Vn|| L2 + ||Vc|| L0O ||Vn|| L2 ||An|| L2 

+ IMIl« II v2c IIl3 + H Vc IIl- IMIl« l|Vc|| L a ||An|| i2 . 

Finally, we show that u G H^Lf* n H 2 L 2 . Testing — Aw to the equations, we have 

^ ||Vo;|| 2 2 + IIV^H 2 .;, < ||-u|| LP \\Vu\\ L i ||V 2 w|| l2 

+ ||Vu|| L 4 ||w||jr,4 ||V 2 a;|| L2 + C \\Vn\\ L2 ||V 2 a;|| L2 , 
where 3 < p < oo and 1/p + l/l = 1/2. Note that, via the Gargliardo-Nirenberg's inequality, 

5 7 

||V«|| L 4 ||w|| L 4 || V 2 w|| i2 < C ||w|| 2 4 ||V 2 w|| L2 < C ||w||| a ||V 2 a;||| 2 . 

We treat the term ||u|| LP HVwH^i ||V 2 w|| L2 similarly to ||u|| iP \\w\\ L i ||Vcj|| l2 in the estimation 
of (12.13j) . Therefore, since V 2 cj G L 2 L 2 , we have 



/ ||Vn|| ioo di< / \\Vu\\ H 2dt<C [ 
Jo Jo Jo 



T 

\uj\\h2 dt < oo. 



This completes the proof. rj 
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3 Global solutions in two dimensions 



In this section, we provide the proof of global existence of smooth solutions in time with large 
initial data in two dimensions. For the proof of Theorem [31 we show some a priori estimates, 
which are uniform until the maximal time of existence. Moreover, such estimates imply that 
the blow-up condition quantity in Theorem [2] is uniformly bounded up to the maximal time 
of existence. Therefore, the maximal time cannot be finite. Now we present the proof of 
Theorem 

Proof of Theorem 3. We first present the following estimates for the solutions to the 
two-dimensional chemotaxis system coupled with the Navier-Stokes equations. 

n(l + |x| + | lnn|) G L°°(0, T; L 1 (IR 2 )), Vy/n G L 2 (0, T; L 2 (M 2 )), (3.1) 

c£ L°°(0, T;L 1 (R 2 )nL 0O (lR 2 )ni7 1 (]R 2 )), Vc G L 2 (0, T; L 2 (IR 2 )), (3.2) 
u G L°°(0, T; L 2 (R 2 )), Vu G L 2 (0, T; L 2 (R 2 )). (3.3) 
We have the mass conservation for n(t,x) as 



n(t, x)dx = I rio(x)dx. (3.4) 

Multiplying c q ~ 1 (t,x) to both sides of the second equation of (jl.ip and integrating over R 2 , 
we have 

-4l|c|lL + 4(g 7 1) ||Vcf|| 2 2 + f k (c)nc"- 1 dx = 0. (3.5) 
q at q J R 2 

Hence, we have c G L°°(0, T; L q ) for any 1 < q < oo and Vc2 G L 2 (0, T; L 2 ) for any 1 < q < oo. 
Multiplying Inn to both sides of the first equation of (jl.ip and integrating over M 2 , we have 



d 
dt 



/ n In ndx + 4 / |V-v/n| 2 (ix + / }c {c)\V c\ 2 ndx = — I x(c)Acn<ix. (3-6) 

JR2 J]R2 J R 2 J R 2 

Multiplying —Ac to both sides of (jl.ip and integrating over M 2 , we obtain 

'll^ 7c ll|2 + II^ c IIl2 = / k(c)Ac ndx + / cd^Uidid^cdx 

7r2 ^ y M 2 



< / A;(c) Ac ndx + (711^^11^2 ||c||loo || Ac||^,2. 

JR2 

Multiplying /x to both sides of (|3,7p and then adding (|3.6p . we have 

/ nlnn + fi\Vc\ 2 dx + / 4| V\/n| 2 + //| Ac| 2 dx 

7r 2 JR2 



(3.7) 



d 
~dt 



< e||Ac|U2||V^Hi4 + C 2 ||c||| 0O ||Vn|| 2 2 + ^||Ac|| 2 2 
< ^||Ac|| 2 2 +eC 3 ||V^|ll2 +C 2 ||c||| 0O ||Vn|| 2 2 , (3.8) 
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where we used the condition (A). Here we choose e to be so small that eC^ < 2 and, for 
convenience, set -4- := C^Hcollf^x,. On the other hand, multiplying u to both sides of the third 
equations of (jl.ip and integrating over K 2 , we have 



Id 

2 dt 1 



H\h + ll Vn lll 2 = - / nVcpudx. (3.9) 
Multiplying (p to both sides of the first equation of (jl.ip and integrating over M 2 , we have 
— / ncj)dx = — I u ■ \7n<f)dx — f Vn • X7(f>dx + J x(c)?iVc • V<f)dx 

dt Jr2 J R 2 J R 2 J R 2 

<-/ u- Vn0dx + C 4 ||VVn|| L 2||Vn|| L 2 +C 5 ||n|| L 2||Vc|| L 2. (3.10) 
Summing (|3.9p and (|3.10p . we have 

sG || * + L^) +||vu| 



2 

II 2 



< CsIIV^IMIVv^IIl 2 +C 6 ||v^|| L 2||Vv^|| L 2||Vc|| l2 . (3.11) 
Multiplying Ai to both sides of ([3. lip and adding (|3.8p . we obtain 

d f A f A 

— / n\nn + ^i\Vc\ 2 + ^-\u\ 2 + \ 1 n<pdx + / 2|V^/n| 2 + ^|Ac| 2 + -±\Vu\ 2 dx 

dt J R 2 2 J R 2 2 2 

< AiCsllnollliUVv^IlL 2 +AiC 6 ||n ||| 1 ||V v ^||L 2 l|Vc|| L 2 

< ^IKIU* + ^IKIM|Vc||| 2 + ||Vv^|| 2 L2 . (3.12) 
Using Gronwall's inequality, we have 



sup ( f n\un + /x| Vc| 2 H -\u\ 2 + Xiruftdx 

o<t<T ViM 2 



+ j T [ |V^| 2 + -|Ac| 2 + — \Vu\ 2 dxdt < C(T). 
Jo Jm? 2 2 

Next, we show that n|lnn| 6 L°°(0, T;L 2 (]R 2 )), following a typical argument for dealing with 
kinetic entropy (see e.g. [5])- We first note that 

/ rc(lnn)_ < C + C I n(x), (3.13) 

2 1 

where (lnn)_ is a negative part of lnx and (x) = (1 + \ x\ ) 2 . Indeed, setting D\ = {x : n(x) < 
e~' x '} and D2 = {x : e - ^' < n{x) < 1}, we have 



/ n(lnn)„ = — nlnn — / n In n < C \/n + / n(x) < C / e '2' + / 



n(x). 
(3.14) 
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n(.,t)lnn(.,t) +/ u||V C (t)||i 2 + ^||^)" 2 



This deduces the estimate (|3.13p . Next, integrating (|3.12p in time t, we get 

\ L 2 \\n(t)(f)\\ L i 

+ f I 2 ( |V^| 2 + ^|Ac| 2 + ^|Vn| 2 ^) dxdr < C 7 + C 8 t + C 9 f \\Vcfr 2 dxdT, (3.15) 
Jo Jr 2 V 2 2 / Jo 

where Cj = J* R2 lnno+// llVcoll^+^y- ||uo||^2+Ai ||no^||xi- Remembering (|3.13j) . we compute 



d 
di 



(x)ndx 



nuV{x)dx+ \ nA{x)dx + / x( c ) n ^cV(x)(i2;. (3.16) 
Vr 2 Jr 2 



The term J" K2 nuV{x)dx is bounded as follows: 



nuV(x)(ix 

Noting that \V (x} \ + \A(x}\ < C, we get 



< IIV^&IMU* < ^V^\\ 2 L 2 +C||n |Ui||u|| 2 2 . 



nA{x)dx 



+ 



x(c)nVcV(x}dx 



< C + C||VV^||r 2 ||Vc|| L2 , 



where we used that ||«|| i2 < C ll^o || ?i II^V^IIl 2 - ^ n summary, we obtain 

4 / {x)ndx < 5\\V^n\\ 2 L 2 + C||u|| 2 2 + C 5 ||Vc||^ + C, (3.17) 

where 5 is sufficiently small, which will be specified later. Therefore, integrating (|3.17p in time, 

(x)n(-,t)dx< f (x)n dx + 5 [ || V>/n|||a + C f \\u\\ 2 L2 + C s [ \\Vc\\ 2 L2 + Ct. (3.18) 
! Jr 2 Jo Jo Jo 



Now adding 2 J n(lnn)_ to both sides of (I3.15p . we obtain 



n(-,t)\\nn(-,t)\+ fi\\Vc{t)\\l 2 + — \\u 



, ^1 ll mil 2 
L 2 



+ Ai lln 



L 1 



+ [ [ (\^Vn\ 2 + fi\Ac\ 2 + \t\Vu\ 2 ) dxdr 
Jo Jr 2 

<C + Ct + C [ \\Vc\\ 2 L2 dxdT + C [ \\uf L2 dxdT, 
Jo Jo 



(3.19) 



where 5 in (|3.17p is so small that term ||V-\/rt|| 2 2 is absorbed to the left hand side of (j3. 15|) . 
Since (I3.19P holds for all t until the maximal time of existence, due to GronwalFs inequality, 
we obtain n|lnn| £ L°°(0, T;L 2 (M 2 )). Moreover, again via the inequality (|3.19p . we deduce 
pT]) - (l33]) . 

We note that from the blow-up criterion in two dimensions in Theorem [2l it suffices to show 
that Vc G £ 2 (0, T;L°°(M 2 )) for global existence of smooth solutions in M 2 . We first consider 
the vorticity equation of velocity fields. Taking curl, we have 



d t u) + (u • V)w - Aco = -V ± n • V0, 
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where V 1 - = (—82, d%). If we multiply u to both sides of the above equation and integrate over 
M 2 , then we have 



-I j /» 

+ \\Vu\\h= / nV4>V L udx<C\\n\\ L 4V A 
1 at j R 2 



Hence, we have 

IMl!°°(0,T;L 2 ) + ll Va; lll 2 (0,T;L 2 ) - ^11 n ll 1 2 (0,T;L 2 ) • 

Since \\n\\ L 2 < \\^n\\ 2 L4 < C\\^/E\\ L 2 ||V^||l 2 , we have u G L°°(0, T; L 2 ) n L 2 (0, T; H l ). Next 
we consider the equation of n. Multiplying n and integrating over M 2 , we have 

1 d f 

olilHli 2 + H Vn lli 2 = / x(c)nVcVndx 

2 at j R 2 

= -\ [ V • ( X (c)Vc)n 2 dx <C f \V 2 c\n 2 dx + C [ \Vc\ 2 n 2 dx 
2 Jr 2 Jr 2 Jr 2 

< C(||V 2 c|| L2 + \\Vc\\ 2 L ,) \\n\\ 2 L , < C||V 2 c|| L2 ||n|| L2 ||Vn|| L2 , 

where we used that x is C 1 an d c G L°°(0, oo;L°°), i.e., x( c ) an d x'( c ) are bounded. Due to 
Young's inequality, we have 

|lhlli 2 + l|Vn||i 2 <C||n||i 2 ||V 2 c||i 2 . 

Therefore, via Gronwall's inequality, we have n G L°°(0, T;L 2 ) n L 2 (0, T;H l ). Multiplying 
A 2 c to both sides of the equation of c and integrating over R 2 , we have 



ld_ 
2~dl 



||Ac|| 2 2 + ||VAc|| 2 2 < ||Vu|| i 4||Vc|| i 4||VAc|| L 2 + ||u|| L c ||V 2 c|| i 2||VAc|| L2 - f k(c)nA< 



cdx. 



We note that the last term above is controlled as follows: 



k{c)nA 2 cdx 



< 



fc'(c)Vc- (VAc)ndx 



+ 



k(c)Vn ■ (VAc)dx 



< C||VAc|| L2 ||n|| L 4||Vc|| L 4 +C||Vn|| L3 ||VAc|| i 2 



<e\\VAc\\ 2 L 2+C\\n\\ 2 L 4Vc\\ 2 Li +C\\Vn\\ 2 L 2. 

Hence, we have 

lj t \\M\h + l|VAc|| 2 2 < C||V«||| 4 ||Vc||| 4 + C7||i i ||£ 0O ||Ac||£ a + C||n|| 2 4 || Vc|| 2 4 + C||Vn||£ 2 . 

< C \\u\\ L2 \\Vu\\ L2 || Vc|| L2 ||V 2 c|| L2 + C || Vu\\ 2 L 2 \\Ac\\ 2 L 2 
+C \\n\\ L2 \\Vn\\ L 2 \\Vc\\ L2 \\V 2 c\\ L2 + C || Vn||£ 2 . 

Gronwall's inequality gives c G L°°(0, T; H 2 )nL 2 (0, T; -ff 3 ), which implies via embedding that 
Vc G L 2 (0, T; L°°). This completes the proof. n 
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4 Global weak solution in three dimensions 

In this section we will show the global existence of the weak solutions for (jl.ip in three dimen- 
sions. We start with notations. ^(M 3 ) is used to indicate the closure of compactly supported 
smooth functions in i/ 1 (]R 3 ) and H~ 1 (R 3 ) means the dual space of -?/q(R 3 ). We also introduce 
the function spaces V(R 3 ), V CT (M 3 ), H(R 3 ) defined as follows: 



3 ) = {u = (u u u 2 ,u 3 ) | m G H^(R 3 )}, V a (R 3 ) = {u£ V(R 3 ) |divu = 0}, 

U(R 3 ) = the closure of V ff (M 3 ) in (L 2 {R 3 )) 3 . 

The dual space of V(R 3 ) is denoted by V'(1R 3 ) = {it = 1*2,1*3) I «i € H~ l (R 3 )}. The 
duality (w,v) for w G V'(R 3 ),v G V(R 3 ) is, as usual, given as (w,v) = Ylf=ii w ii v i)H^ 1 xH^ 
and we denote V°(R 3 ) = {w G V'(R 3 ) | (w,v) = for all v G V CT (M 3 )}. 
Next, we define the notion of a weak solution for the system 

Definition 5 Let < T < 00. A triple (n, c, it) is called a weak solution to the Cauchy 
problem (jl.ip in R 3 x [0, T) if the following conditions are satisfied: 

(a) The functions n and c are non-negative and (n, c, it) satisfy 

n(l + \x\ +|lnn|) G L°°(0, T;L 1 (R 3 )), Vy/n G L 2 (0, T;L 2 (R 3 )), 

c G L°°(0, T; L^IR 3 ) n L°°(M 3 ) n ^(M 3 )), Vc G L 2 (0, T; L 2 (R 3 )), 
it G L°°(0, T; L 2 (M 3 )), Viz G L 2 (0, T; L 2 (R 3 )). 

(b) The functions n,c, and it solve the chemotaxis- fluid equations (II. ip in the sense of dis- 
tributions, namely for any * G C^QO, T]; (C^°(M 3 )) 3 ) untfi V • * = 

I (it •*)(•, T) + [ [ u ■ (dt& + A*) + f [ i*®u:W 
JR 3 Jo JR 3 Jo JR 3 

poo r r 

- / nX7(p-^+ «o-*(0,x) = 0, 

Jo JR 3 JR 3 

where u ® it : = ^ fc=i u^u k dj^ k and 

/ / n(d t cp + Acp) + / nu-V(p+ / x(c)nVc-Vip+ n (x)(p(0, x) = 0, 
JO Jr 3 JO Jr 3 JO Jr 3 Jr 3 



c (xM0,x) = 



/ / c(a^ + Av?)+ / cu-Vcp- / fc(c)ny»+ / 

Jo Jr 3 Jo Jr 3 Jo Jr 3 Jr j 

for any ip G C^T]; (C^°(1R 3 ))) with <p(;T) = 0. 
(cj T/ie functions n, c and u satisfy the following energy inequality: 
r U|2 I Vcl 2 /" T 

/ ( i ^-+n0 + n|lnl^| + ^-^ + (x)n)c^x+ / \\Vu\\ 2 L2 + \\VVn\\ 2 L 2 + \\Ac\\ 2 L2 dt < C, 
Jr 3 * ^ Jo 

with C = C{T, ||x(c)||loo, ||(x)n ||Li, ||Vc ||z,2, ||n | lnn |||ii, ||A0||i<x., ||V0||l°°, ||0||z,°°). 
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Now we compute a priori estimate of an energy inequality under the Assumption (AA) 
and (B). We note first, by maximum principle, that 

n(t,x)>Q, c(t,x)>0, \\c(t)\\ LP < \\co\\lp for t > 0, 1 < p < oo. 

It is straightforward that = H^oIIl 1 f° r t — an d 

d ( f In' 2 



-dx + / ri(j)dx ) + 

dt \./m>3 2 



\X7 u \ 2 dx= / nAcpdx+ / x(c)nV cV 4>dx , (4.1) 

— / n In ndx + [ Lj x _|_ /" ^'(c)|Vc| 2 ncix = — f x(c)Acncfa;, (4-2) 

Jr 3 Jr 3 n Jr 3 Jr 3 

\Vc\ 2 dx + / |Ac| 2 dx = / k(c)Acndx + / cdidjcdiUjdx. (4.3) 

JR 3 y R 3 ^ J R 3 



d 

dt 

Multiplying u to the last equation fj4.3|) and adding it to the second equation ()4.2p . we have 

^- (/ nlnndx + fi\Vc\ 2 dx] + [ \V\/n\ 2 dx + /i / | Ac| 2 dx + / x'(c)| Vc| 2 ndx 
a* Vjr 3 / JR 3 JR 3 J 

f c 

<- / (x(c)-M*(c))Acndx + M||co||L<»||Vu|| £a ||A C || La < ^||Vn|| 2 2 + ^||Ac|| 2 2 (4.4) 

Jr 3 v > ' 2 4 

=0 

for some C\, which can be taken bigger than 1, i.e. C\ > 1. Also it holds that 

/ {x)ndx = [ nuV(x)dx + f nA(x)dx + f x( c ) n ^cV '(x)dx. (4.5) 
dt Jm. 3 Jr 3 Jr3 Jk.3 

Since the term J K3 nuV{x)dx is bounded as follows: 

\\n\\ L e\\u\\ Lfi < Cllnll^UVv^llillVnll^ < ^||Vv^||£ 2 + ^||V«||| 2 + C(||n || L i), 
we can have 

/ (x)ndx< i||Vv^|| 2 2 + i||Vn|| 2 2 + / X (c)nVcV(x)dx + C. (4.6) 
at j R 3 I l j R 3 

We estimate the term _L 3 x(c)nVcV(x)dx similarly as above. 

X (c)nVcV(x)dx < CUv^HiallVclUa < C\\M\l* II Vy^U* II Vc||| 2 1| Ac||| 2 

3 

< C||V C || L2 ||Ac|| L2 + iyVy^llia < C||Vc||| 2 + i||A C || 2 2 + i||Vv^||i a . ( 4 -7) 



Multiplying C\ to (14. ip and adding it together with (|4.4p and ()4.6p . we have 

/ /" |u| 2 |Vc| 2 N 

I / - — h ri(j)) + n In n -\ \- (x)ndx 

\Jr 3 ' 2 2 



r/ 

+ ^^i||V«|& + J||Vv^|| 2 i2 + i||Ac||| 2 < C(||Vc|| 2 2 + \\uf L2 ) + C. (4.8) 
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Then, by GronwalFs inequality, we have 

12 |V7„|2 



/ ( 1- n<p + n in ra H \- \x)n)dx 

Jr 3 2 2 

+ / l|Vn||2 2 + HWiTHia + ||Ac||| 2 dt < C, (4.9) 

where C(T, ||x(c)||x,°o, ||no||z,i, ||(x)no||z,i, || A</>||x,°o, ||V0||ioo). By same reasoning for treating 
n(lnn)_ term in (|3.14p . it follows that 

[ (\u\ 2 +ri(j) + n\lnn\ + \Vc\ 2 + {x)n)dx + [ \\Vu\\\ 2 + y/n\\\ 2 + \\Ac\\\ 2 dt < C. (4.10) 
JR 3 Jo 

Streamline of constructing global weak solutions, as in usual steps for the Navier-Stokes equa- 
tions, is the following: 

■ regularizing the system for which we prove the existence of smooth solutions 

■ finding uniform estimates for the solutions of the regularized system 
• passing to the limit on the regularized parameters. 



4.1 Regularization 

In this subsection, we intend to construct approximate solutions of the system. For the 
incompressible Navier-Stokes equations defined on a general bounded domain, the global weak 
solutions are constructed by using the spectral projections (Pk)kez, associated to the inhomo- 
geneous Stokes operator ([TJ Chapter 2]). A number of useful properties of the family (-Pfc)fcez 
are listed as follows: For any u £ T-L{fl), 

PkPk'U = P min (k,k>) u , , lim \\ p kU - u\\ H(n) = 0, (4.11) 

||VP fc u|| L 2 (Q) < v^llull^), \\AP k u\\ L 2 (n) < &|H| L 2 (n) , (4.12) 

\\{l-P k )u\\ L 2 < ^=Nk- (4-13) 

In particular, (|4.12p implies P k u G L°°(S7) for u G L 2 (0) in three dimensions. 
Definition 6 The bilinear map Q is defined by 

Q : V x V -> V, 

(u, v) i — y — div (u ® v). 

From now on we denote by Ti k (M?) the space -Pfc%(R 3 ). We regularize (jl.ip by a frequency 
cut-off operator Pk and a mollifier a e : 

' d t n k ' e (t) = -u k ' e ■ Vn k ' e + An k ' e - V • (n k > e [(x(c k > e )Vc k > e ) * a e }), 
< d t c k ' e (t) = -u k ' e ■ Vc k ' e + Ac k ' e - k(c k ' e )(n k ' e * a e ), (4.14) 
k d t u k > € (t) = -P k Q(u k ' e , u k > e ) + P k Au k > e - P k (n k ^V<p), 

with initial data 

/ k,€ k,e fc,e\ / e en e\ 

(n ' , c ' , u ' ) = (n * a , c * a , P fe u * (J ), 
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where uq,cq,uq is the initial data of (jl.ip satisfying the condition (jl.8p in Theorem [H The 
mollifier is defined as usual such that a e (x) = e~ s a(e~ 1 x) for a £ Cq°(R 3 ). Apart from the 
frequency cut-off the regularization is same one for a chemotaxis-fluid model studied in [13| . 
Repeating similar arguments in Theorem [TJ we obtain the local solution of (jl.ip in the class 

c k,e e L °o( 0) T; J H' ?n - 1 (]R 3 ))nL 2 (0,T; J ff m (]R 3 )) (4.15) 
G L oo(0 ) r; J ff m - 1 (R 3 )nHifc(M 3 ))nL 2 (0,r;iJ m (R 3 )) 

for some time T and for all m > 3. It turns out that due to the regularization of nonlinear 
terms and smoothing properties of Pk (see (I4.12p ). the local solution of (11.11) can be extended 
up to infinite time. 

Proposition 2 The regularized system (|4.14p has the unique global solution (n k,e ,c k,e ,u k ' e ) in 
a class (|4.15p for any time T < oo. 

Before presenting the proof we observe that the approximating solution (n k,e , c k,e , u k ' e ) of (|4.14p 
satisfies an energy inequality. 

Proposition 3 The solution (n k ' € ,c k,€ ,u k,e ) of (14.14p satisfies the following inequality. 

,k,e\2 |V7„fe,e|2 



\ (- + n k ^4>) + n fc ' e | lnn fc,e | + + (x)n k ' e dx 

Jr3 2 2 

rT 

+ / ||Vu fc ' e ||| 2 + ||VVn^||| 2 + \\Ac k ' e f L2 dt < C, 
Jo 

where C = C (T, \\x(c)\\l°° ,\\(x}no\\ L i,\\Vco\\ L 2, IKI lnno||| L i, ||A0||z,oo, ||V<^|| 



(4.16) 



Proof. We note that the same cancellation as in (|4.4p holds for the regularized system (|4.14p . 
hence (n k,e ,c k,e ,u k,€ ) satisfying (I4.15P satisfy the energy inequalities (|4,ip - (l4.5p . Moreover the 
following moment bound holds by similar estimates as (|4.6p . (|4.7|) . 



d_ 

dJ 



(x)n k > e dx = [ n k > e u k ' e V{x)dx + f n k ' £ A{x)dx+ [ n fc ' e [(x(c fc ' e )Vc fc ' e ) * a e ]V(x)dx 

■ JM.3 J R 3 J R 3 

< C\\u k ^\\ L2 \\Vu k >*\\ L 2+C\\Vc k > € \\ L 2\\Ac k > e \\ L 2 + ^\\VV^M\\ 2 L 2 + \\n Q \\ L i. 



Then we have (I4.16P with T depending on ||VcQ ,<: || i 2, || (a;)n ' e ||i,i, H^o'l ^"o'lllx 1 - ^ ^ s 
immediate to have 

II Vc o'1l 2 + IK^WIU 1 ^ II Vc oIIl2 + 11(^)^01^1- 

Note that xlnx is convex and dfi = a e (y)dy provide a probability measure. Then by Jensen's 
inequality, we have 

n^ e (lnn^ € ) + < (n (lnn )+) * a € . 

Integrating the above in x and observing that lim^o ||(no| lnno|) * o" e || L i = ||no| hmo| H^i, we 
have 

||no' e (lnno' e ) + || L i < ||n | lnn | || L i. (4.17) 



22 



For the ||n ' e (lnnQ' e )_||^i, proceeding similarly as (|3.14p . we have 

||4' e (lnng' e )_|| L i <C+ [ n k ' t (x)dx < C ( 1 + / n(x)dx) , 

JR3 \ Jm 3 J 

from which we deduce the proposition. rj 

Now we give the proof of Proposition [2j 

Proof of Proposition [2] We first observe that the regularity criterion in Theorem [2] hold 
true for the system (|4.14p . Since its verification is tedious repetition of that of Theorem [21 
we omit its details. If we consider the second equation of (I4.14h . then we have the following 
energy estimates. 



1 d \\r k ' e \\ 2 -4- \\X7r k < e \\ 2 < r\\X7n k > e \\ 2 \\r k > e \\ 2 A- CWX7 r k,e \\ 2 \\n k ^\\ 2 4- 1 II X7r k ' e II 2 

2~di<< Hh 2 + H Vc \\h 2 - u W Vu Hl 3 H c Hl 6 + L/ H Vc llL 6 ll n Mz, 1 2 " " H 



< C^ll^ll.llc^ll^ +C||n^||| 1 ||c fe -l| 2 + ^||Vc fe -l| 2 . 

By using GronwalPs inequality, we have ||Vc fc,<E || roor2 < 00. Since ||w fc ' e (i)||z,2 is bounded and 

||Vn fc,<E || i 2 < Cv^II^^^IIl 2 i we can a ls° demonstrate that the Serrin condition in Theorem [2] is 
satisfied for u k ' e . This completes the proof. □ 

4.2 Global weak solutions 

In this subsection, we give the proof of Theorem [U 

Proof of Theorem 4. We consider an approximating sequence (n ' e , c ' e , Uq € ) to (iio,cq,Uq). 
Note that 

/ \tiq — no I + |Vc ' e — Vco\dx + / |u ' e — u\ 2 dx — > 0., 

7R3 Jm.3 



and 



/ (x)n l Q e dx+ / n ' e | In n ' e | < C \ (x)nodx + / no| \nno\dx + C. 
Jm.3 Jm.3 J K 3 j M 3 



We denote by (n l,€ , c l,e ,u l,e ) the approximating solution constructed in the previous section for 
the system (|3~T4l with initial data (W(0, ■), c'(0, ■)) = (n (0, •), c (0, •)) and u'(0, •) = fjito(-)- 
Several uniform estimates hold for the approximating solutions: 

II c/ '1Il°°(o,T;Lp(r3)) < C for 1 < p < 00, (4.18) 



|L°°(0,T;LP 

II^^IIl^^T;^ 1 ^ 3 )) + II^ c/ ' £ |Il 2 (0,T;L 2 (R3)) < C, (419) 

^Hl^t^CrS)) + l|V\ / n 7 ^|| L 2 (0iT;L 2 (IR 3 )) < C, (4.20) 



e "L°°(0,T;L 2 (R 3 )) + II ||l 2 (0,T;L 2 (R 3 )) < C- ( 4 -21) 



\U 



Then there exists subsequences n ,e , cr ,e ,u > e and some functions n,c,u such that 



n 1 



^^v 7 ™ L°°(0,r;L 2 (M 3 )) - weak*, 



c'' e c L°°(0,T;L p (M 3 )) n L°°(0, T; tf^lR 3 )) - weak*, 
u i,e _± u L oo^ 0) T . L 2( M 3^ _ weak * n L 2( 0j T . V CT (M 3 )) - weak 
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for 1 < p < oo. Let us show that n, c, u is a weak solution in the sense of Definition By 
Gagliardo-Nirenberg inequality and (|4.20|) . we have 

\n l ' £ \ p dx < C||n || i ? (]R 3)||Vv / n^||^ 2 p (IK 3 ) ) , 



and therefore, 

II^'1l 9 (o,t ; lp(r3)) < C(T), 1 < q < (4.22) 

for 1 < p < 3. Some strong convergences are necessary. We note that (|4.22p implies the source 
term of the Navier Stokes equation n l,€ 'V<p is in L 2 ([0,T];V a '(M. 3 )) uniformly with respect to I; 
for any w £ L 2 ([0, T]; V cr (lR 3 )), it holds that 



T 

/ Pi(n l ' e V<p)wdxdt < IIV^II^^)!!^' 6 !^ otl | ||w|| L 2 (0iT;L 6 

It proves that dtu l ' e is uniformly bounded in L 2 (0, T; V</(M 3 )). Note that is uniformly 
bounded in L°°(0, T;Ti(R 3 ))) n L 2 (0, T; V CT (]R 3 )) due to (ITOjl . Combining these facts and 
(|4.12p . (|4.13|) we have compactness result for (u l,€ ) (see [TJ Proposition 2.7] for detailed proof): 
there exists u in L 2 (0, T; V<j(K 3 )) such that up to subsequence 

lim / / \u l ' e (t,x) - u(t,x)\ 2 dxdt = 0, (4.23) 

l-x>o,e-+0j J K 

for any T > and compact subset K of M 3 . In addition, for ^ G L 2 ([0, T]; V(M 3 )) and 
$ e £ 2 ([0,T] x M 3 ) 



lim / / Vu l > e (t,x)V$(t,x)dxdt = [ [ Vu(t,x)V^(t,x)dxdt, 

/^oo,£^0j J R 3 J J R 3 



JR J J 

lim / / u l,t (t,x)$(t,x)dxdt = / / u(t,x)§(t,x)dxdt. 



(4.24) 



Furthermore, For any ^ £ C 1 (M + ;V 



lim sup 

Z-^oo,e->0 ie [ 0]T ] 



(u l,e (t, x) — u(t, x))ip(t, x)dx 



0. (4.25) 



(4.26) 



Applying a test function ^ in C^QO, T]; V CT (M 3 )), we obtain 

^(u l > e (t)Mt)) = <A^(t), W)> + (Q(u^(t),u^(t)),m(t)) 
+ ((n'' e V^),fl*(t)> + (u , ' e (t),^*(t)>. 
Following the arguments in [I], that is, using ()4.23p - ()4.25p and the fact 

lim sup ||fi*(*)-*(t)|| v(R s ) = 0, (4.27) 
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we can pass to the limit with respect to I so that 

/ u ■ ^(T, x)dx + ( j (Vn : - u<g> u :W - u ■ d t ^)(s,x)dxds 
Jr 3 Jo Jr 3 

u (x)^(0,x)dx + lim / {n l ' e V(j),^)dt. 

For the strong convergence of (n l,e ) we have VttM — > ^/n strongly in Lf oc (M. 3 ) for a.e. t G [0, T] 
by Sobolev embedding. Since ||Vn Z ' e (i)|| L 2( R 3) is continuous in time, we redefine n(t) such that 

||Vn'' e — v / ^ll-L 2 (K 3 ) ~~ * f° r an ^ ^ [0, 2^]. Then by (|4.22p and Lebesgue Dominated convergence 
theorem, it follows that 

; 2p 

II"'' 6 " n\\ Lq{0 ,T;Ll c m) ^0, 1 < q < (4.28) 

for 1 < p < 2. For the convergence of (c'' e ) we have c l,€ (t) — > c(t) strongly in Lf oc (M. 3 ) for all 
t G [0, T] and therefore, 

H c/ ' e - c llif oc ((o,T)xR3) ->0, l<p<oo (4.29) 
by the uniform boundedness (|4.18p . Moreover we have 

||Vc'< e - Vc|| L2(0iT . LL(R3)) 0, 1 < p < 6. (4.30) 
By (I4.19p . ||Vc'' e || L 2( 0)T; ^i( R 3)) is uniformly bounded. For any Vg G L 4 (0, T; L 2 (M 3 )), we have 

/ / d t Vc l ^gdxdt< [ [ u l ' e Vc l ^Vg + Ac l > e Vg + k(c l > e )(n l > e *a e )Vgdxdt. 
Jo Jr 3 Jo Jr 3 

We estimate 

/ u l ' e Vc l > e Vgdxdt<C i ||n i 'l L6(R 3 ) ||V C i 'l| 2(R3 JAc z ' £ ||| 2(R 3jV(7|| L 2 {R 3 )( it 
o Jr. 3 Jo v ' v ' 



f T - 



< C||Vn i ' e || i 2( 0)T . i 2( R 3))||Ac'' e ||2 2 ( 0iT; L2( R 3))||V5(|| L 4( 0iT;£ ,2( R 3)), 

lo L k( - cl ' £ ^ nl ' £ * a ^ V 9dxdt < C II (0 ^.^ (R3)) II Vffll i4 (0 ^ ;jL2(R3)) . 

Thus we have d t c l ' e G L5 (0, T; #~ 1 (M 3 )). The strong convergences (I4.28p - (|4.30p are enough to 
pass to the limit for nonlinear terms in the chemotaxis part. For instance, testing a Vl/ E C^°(IR 3 ) 
to the worst nonlinear term V • (n l,e (x(c l,£ )Vc l ' e ) * a e ), we have 

/ / V ■ {n l > e [(x(c l ' e )Vc l ' e )*a e })^ - V ■ (n X (c)Vc)^dxdt 
Jo Jr 3 

T r 

(n l ' e - n)[(x(c /,£ )Vc'' e ) * a t ]V^dxdt 



i 



o 



+ 



/ / n[(x(c l ' e )Vc l > e ) *a 6 - X (c)Vc]V^dxdt. 
Jo Jr 3 
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The second integral is 

/ / [(nVf) *a e - n^]x(c /,£ )Vc'' e + nW(x(c'' e )Vc'' e - X (c)Vc)dxdt 
Jo Jk.3 

= [ [ [(nV^)*a e -n^}x(c l ^)Vc l ^ + nV^([(x(c l ^)-x(c)]Vc l ' £ + x(c)(Vc l ^-Vc))dxdt. 
Jo Jm. 3 

The integrals go to zero by the uniform estimates (|4.18|) - (|4.2ip and (|4.28p - (|4.30|) with the Lip- 
schitz continuous assumption on x(")- 

Lastly, we consider the approximated energy inequality (|4.16p replacing n l,e \ lnn^l with n l ' £ lnn'' e . 
Taking the limit and uing the convexity of x In x we deduce 

r | M |2 iy c |2 pT 

/ (-^f" +n<f) + nlnn+ + (x)n)dx + / \\Vu\\\ 2 + \\V ^\\\ 2 + \\Ac\\\ 2 dt < C, 

Jm 3 2 2 Jo 

with C = C(T, ||x(c)||l°°, ||(g)ra || L i, ||Vcq|| £ 2, ||n | lnn ||| L i, ||A0||x,oo, ||V0||j,oo, ||<^||l°°)- By 
the same reasoning for treating ?i(lnn)_ term in (|3.14p we show the weak solutions (n, c, u) 
satisfy the energy inequality in Definition [5] (c). This completes the proof of Theorem 4. m 

Acknowledgments 

M. Chae's work was supported by the National Research Foundation of Korea(NRF No. 2009- 
0069501). K. Kang's work was partially supported by KRF-2008-331-C00024 and NRF-2009- 
0088692. J. Lee's work was partially supported by NRF-2009-0072320. We appreciate professor 
Dmitry Vorotnikov for valuable comments. 

References 

[1] J.Y. Chemin, B. Desjardins, I. Gallagher, and E. Grenier, Mathematical Geophysics, Ox- 
ford university press, 2006. 

[2] A. Chertock, K. Fellner, A. Kurganov, A. Lorz, and P. A. Markowich, Sinking, merging 
and stationary plumes in a coupled chemotaxis- fluid model: a high-resolution numerical 
approach, 2011, to appear in J. Fluid Mech. 

[3] L. Corrias, B. Perthame, H. Zaag, A chemotaxis model motivated by angiogenesis, C. R. 
Math. Acad.Sci. Paris, 336(2); 141-146, 2003. 

[4] L. Corrias, B. Perthame, H. Zaag, Global solutions of some chemotaxis and angiogenesis 
systems in high space dimensions, Milan J. Math., 72(1); 1-28, 2004. 

[5] R. Diperna and P. L. Lions, On the Cauchy problem for Blotzmann equations: Global 
existence and weak stability, Ann. Math., 139: 321-366, 1989. 

[6] R. Duan, A. Lorz, and P. Markowich, Global solutions to the coupled chemotaxis-fluid 
equations, Comm. Partial Diff. Equations, 35(9); 1635-1673, 2010. 

[7] M.D. Francesco, A. Lorz, and P. Markowich, Chemotaxis-fluid coupled model for swim- 
ming bacteria with nonlinear diffusion: global existence and asymptotic behavior, Discrete. 
Conti. Dyna. Sys., 28(4); 1437-53, 2010. 



26 



[8] Y. Giga and H. Sohr, Abstract LP estimates for the Cauchy problem with applications fo 
Navier-Stokes equations in exterior domains, J. Funct. Analysis, 102(1): 72-94, 1991. 

[9] M.A. Herrero and J.L.L. Velazquez, A blow-up mechanism for chemotaxis model, Ann. Sc. 
Norm. Super. Pisa, 24(4); 633-683, 1997. 

[10] D. Horstman and G. Wang, Blow-up in a chemotaxis model without symmetry assump- 
tions, European J. Appl. Math. 12(2); 159-177, 2001. 

[11] E.F. Keller and L.A. Segel, Initiation of slide mold aggregation viewd as an instability, J. 
Theor. Biol., 26(3); 399-415, 1970. 

[12] E.F. Keller and L.A. Segel, Model for chemotaxis, J. Theor. Biol., 30(2); 225-234, 1971. 

[13] J.-G. Liu and A. Lorz, a coupled chemotaxis-fluid model, Ann. I. H. Poincare, 28(5): 
643-652, 2011. 

[14] A. Lorz, Coupled chemotaxis fluid model, Math. Models and Meth. in Appl. Sci., 20(6) :987- 
1004, 2010. 

[15] T. Nagai, T. Senba, and K. Yoshida, Applications of the Trudinger-Moser inequality to a 
parabolic system of chemotaxis, Funkcial Ekvac. 40(3); 411-433, 1997. 

[16] K. Osaki and A. Yagi, Finite dimensional attractors for one- dimensional Keller-Segel 
equations, Funkcial Ekvac. 44(3); 441-469, 2001. 

[17] C.S. Patlak, Random walk with persistence and external bias, Bull. Math. Biol. Biophys. 
15; 311-338, 1953. 

[18] L. Tartar, Topics in nonlinear analysis , Publicatons mathematiques de l'Universite de 
Paris-Sud(Orsay), 1978. 

[19] Y. Tao, Boundedness in a chemotaxis model with oxygen consumption by bacteria, J. Math. 
Anal. Appl. 381(2); 521-529, 2011. 

[20] I. Tuval, L. Cisneros, C. Dombrowski, C. W. Wolgemuth, J. O. Kessler, and R. E. Gold- 
stein, Bacterial swimming and oxygen transport near contact lines, PNAS, 102(7): 2277- 
2282, 2005. 

[21] M. Winkler, Aggregation vs. global diffusive behavior in the higher- dimensional Keller- 
Segel model, J. Differential Equations, 248(12); 2889-2995, 2010. 



Myeongju Chae 
Department of Applied Mathematics 
Hankyong National University 
Ansung, Republic of Korea 
mchae@hknu . ac . kr 



Kyungkeun Kang 
Department of Mathematics 
Yonsei University 
Seoul, Republic of Korea 
kkang@yonsei.ac.kr 



Jihoon Lee 
Department of Mathematics 
Sungkyunkwan University 
Suwon, Republic of Korea 
jihoonlee@skku.edu 



27 



